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Abstract 

Using the fact that pseudospin is an approximate symmetry of the Dirac 
Hamiltonian with realistic scalar and vector mean fields, we derive the wave 
functions of the pseudospin partners of eigenstates of a realistic Dirac Hamil- 
tonian and compare these wave functions with the wave functions of the Dirac 
eigenstates. 
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Pseudospin doublets were introduced more than thirty years ago into nuclear physics 
to accommodate an observed near degeneracy of certain normal-parity shell-model orbitals 
with non-relativistic quantum numbers (n r , £, j = £ + 1/2) and (n r — 1, £ + 2, j = £ + 3/2) 
where n r , £, and j are the single-nucleon radial, orbital, and total angular momentum 
quantum numbers, respectively [0,0. The doublet structure, j = £ ± s, is expressed in 
terms of a "pseudo" orbital angular momentum £ = £ + 1 coupled to a "pseudo" spin, s = 
1/2. This pseudospin "symmetry" has been used to explain features of deformed nuclei ||, 
including superdeformation J| and identical bands |5||J. Although the observed reduction 
in pseudo spin-orbit splitting follows from nuclear relativistic mean-fields 0, only recently 
has the pseudospin "symmetry" been shown to arise from a relativistic symmetry of the 
Dirac Hamiltonian |8|U. 

The Dirac Hamiltonian, H, with an external scalar, Vs, and vector, Vy, potentials is 
invariant under a SU(2) algebra for Vs = Vy + constant, leading to pseudospin symmetry 
in nuclei 0. The pseudospin generators, S^, which satisfy [H , S^} = in the symmetry 
limit, are given by 
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where = are the usual spin generators, the Pauli matrices, and U p = (T p — is 
the momentum-helicity unitary operator introduced in R]. If, in addition, the potentials are 
spherically symmetric, Vsy(r) = Vsy(r), the Dirac Hamiltonian has an additional invariant 
SU(2) algebra, [H , L M ] =0, with the pseudo-orbital angular momentum operators given by 
= ^ 1° ) . Here £ p = U p 4 4 = r x p, while ^ = 1^+4 = ( 4+ 5 m ) = 4+^- 
The eigenfunctions of the Dirac Hamiltonian are also eigenf unctions of the Casimir operator 
of this algebra, L • L \h r , £, j, m) = £{£ + l)|n r , £, j, m), where we have used a coupled 
basis, j = £ + s, and set % = c — 1. Here j is the eigenvalue of the total angular momentum 



operator J M = L M + Sp, J • J \h r , £, j, m) = j(i + 1) |^r, ^, j, m), m is the eigenvalue of 
and n r is the pseudoradial quantum number which we define below. 

In the pseudospin symmetry limit, the eigenstates of the Dirac Hamiltonian in the doublet 
j — £ ± 1/2 are degenerate, and are connected by the pseudospin generators S^. 

Sp\fi r , 1, ji, rrii) = J2 A j s ,m s ,j^mAn r , £, j f , m f ) . (2) 

Here A jf>mfiji>mi = (-i)|-m/4i ^+i)W/+D i where the symbols are 
Wigner 3-j and 6-j symbols respectively However, in the exact pseudospin limit, = —Vy, 
there are no bound Dirac valence states. For nuclei to exist the pseudospin symmetry must 
therefore be broken. Nevertheless, realistic mean fields involve an attractive scalar potential 
and a repulsive vector potential of nearly equal magnitudes, V$ ~ — Vy, and calculations 
in a variety of nuclei confirm the existence of an approximate pseudospin symmetry in 
the energy spectra [JlQ| — p72 ] . Since pseudospin symmetry is broken, the pseudospin partner 
produced by the raising and lowering operators acting on an eigenstate will not necessarily 
be an eigenstate. The question is how different is the pseudospin partner from the eigenstate 
with the same quantum numbers? As noted, energy splittings suggest that the breaking of 
pseudospin symmetry is small, but is the breaking in the eigenfunctions small as well? While 
previous studies have compared the lower components of the Dirac wave functions of the two 



states in the doublet fll0|-|12[1, it is the behaviour of the upper components which is of most 
interest since they dominate the Dirac eigenstates. The relativistic pseudospin symmetry 
has unique and interesting features in the following sense. First, the pseudospin generators 
of Eq. ([I]) intertwine space and spin, and thus lead to an uncommon symmetry structure of 
doublets with different radial wave functions. Second, since bound Dirac valence states do 
not exist in the symmetry limit, the pseudospin properties of realistic wave functions can 
not be determined by perturbation theory. These aspects motivate the present study. 

To determine the pseudospin partners we need to expand the Dirac eigenfunction into a 



spherical basis, 

( r | h r , I j = £ + ~ m) = (g^ hl3 {r)\Y^\r)x] { £ , «/^ i (r)[K®(/ i )xlg>) (3a) 
(r\n r , I 5 = i~\, ™> = (^yMF^ (3b) 

where Y^}(f) is the spherical harmonic and \ is the s pi n function. From this expansion 
we see that the pseudoradial quantum number, n r , is the radial quantum number of the 



lower component of the Dirac eigenfunction [13] as well as the radial quantum number of 



the upper component of the eigenstate with j — £ — 1/2. Because the pseudospin generators 
depend on the unit momentum vector p, we transform the eigenf unctions to momentum 
space in order to calculate the effect of the pseudospin generators on the eigenfunctions: 



P I n r , I j = 1+ i m) = (^ r _ li j )i (p)^ (l+1) (p)x]^^i/^ ) ^-(p)^ (l) (p)x]^ ) ) (4a) 
P I n r , £, 3 = i-\,m) = (^,^>)[r^(p)x]ff ,</« P ,j; i (p)[i r ®(p)x]S?) • (4b) 
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The corresponding spherical Bessel transforms of the radial wave functions are given by 

^ r°° . , s , x o . ~„ l 



/ 2 r°° ~ i 

9n r -l,I,j(P) = (- l Y +1 ^- J 3 i+ M) 9n r -lJ,j( r ) r2dr J=i+ 2 (5a) 

/2 r°° 1 
- J Q Ji^(P r ) 9n r Jj(r) r 2 dr 3 =£-- (5b) 

/,,,.,(/') = (-iY\j- I 3-M) fnrjjr) r'dr j = £ ± -. (5c) 

We then are able to derive 

^ fi\ Tiri £■> jii 771^ ) ^4j'j,mi,ji,mi | 7l ri £■> jii 7Tl{ ) ~\~ ■ / ^-jf,mf,ji,nii I ^r? ^> jfi TTlf )^ ^ . (6) 

Here the superscript psp on the second term denotes the pseudospin partner with jf ^ ji. 
Even with pseudospin breaking, the pseudospin generators do not change £. In addition, 
from Eq. @ we see that the first term with jf = ji is exactly equal to the original eigenstate, 
independent of the amount of pseudospin breaking. This follows from the orthogonality of 



the spherical Bessel functions, ^ J °° p 2 dp j-{pr)j-{P x ) — ^r^- F° r the partner with jj ^ j i: 
the wave function in coordinate space reads 

(r \n r , I, j = t+\, m)^= (^_ ll^l /r)[r( ^+1 )(Ox]S^ i/^/O^COx]^) (7a) 
(r\n r , I j = i-\, m)*"= (^-/r)[r( ^ - 1 )(Ox]S^ i/n.iW^^x]^) ( 7b ) 

where in general n' r can differ from n r since the states with jf ^ ji in Eq. (|6|) are not 
Dirac eigenstates. The expressions for g psp ~ t (r) with j = £ + \ and g psp ~ (r) with 
j=i-\ involve a double integral - J J™p 2 dpJ™ x 2 dx over j l _ x {pr)j- +x {px)g hr _ l j j {x) and 
i i+1 {pr)j i _ i {px)g flr i ■ (x) respectively. The p-integration can be carried out and altogether 



we find 14 



C-Ui+iW = 9, r ,ijr) - ^1 [ dx x^g^x) j=l-\ (8b) 

In the pseudospin limit (Vg + Vy= constant) 

|ra r , £, j, m) psp = | TV, I, j, m ) . (9) 

Since pseudospin symmetry is slightly broken in nuclei, the pseudospin partner can differ 
from the Dirac eigenstate and it is of interest to examine the deviations from the condition of 
Eq. @. Dirac bound states satisfy g hr ~ t j = i± 1 / 2 ~ r£±1 f° r small r, and fall off exponentially 



~ exp(— \/ M 2 — E 2 r), for large r [13|, where M is the nucleon mass and E is the Dirac 
eigenenergy. Consequently, as seen from Eq. (8a), for the Dirac eigenstate with j = £ + 1/2, 
its pseudospin partner g p f p ^ J _ 1 has the expected behavior for small and large r. On the 
other hand, as seen from Eq. (8b), for the Dirac eigenstate with j = I — 1/2, its pseudospin 
partner g psp _ . ~ r £_1 for small r, and falls off as a power law r~^ +2 ^ for large r. As such 
it has a behavior which is very different from that of a Dirac bound state with j = I + 1/2. 



This asymmetry in the behavior of the pseudospin partners of j = £+1/2 or j = £ — 1/2 
Dirac eigenstates, is evident in the analysis of nuclear wave functions presented below. These 
realistic wave functions were obtained in a relativistic point coupling model, and we refer 
the reader to [K| for details on the parameterization of the potentials, and the data that 
has been used to fix it. 

We first examine Dirac eigenstates with j = i + 1/2 and wave functions as in Eq. (3a). 
Their partners with j' = j — 1 are obtained from Eqs. (8a, c). As an example, we consider 
the realistic relativistic mean field Dirac eigenstates 0d 3 /2, 1^3/2 — 1, j — 3/2) for 
208 Pb jKj. In Fig. 1 we compare the spatial wave functions of these pseudospin partners, 
[P(0d 3 /2)]si/2, [P(ld 3 /2)]si/2 with the eigenstates, lsx/ 2 , 2si/ 2 . (The symbol P means the 
S1/2 partner of the 0d 3 /2 or Ig^/2 eigenstates). The lower components agree very well, which 
was noted previously [ji~0|-|l2"|, except for some disagreement on the surface. For the upper 
components the agreement is not as good in the magnitude but the shapes agree well, with 
the number of radial nodes being the same [n' r = h r in Eq. (7a)]. The agreement improves 
as the radial quantum number increases, which is consistent with the observed decrease in 



the energy splitting between the doublets P,[10|- As another example in the same category 
(j = t + 1/2), we compare in Fig. 2 the [P(0h 9 / 2 )]f7/2 partner of the 0/i 9 / 2 eigenstate 
(£ = 4, j = 9/2) with the I/7/2 eigenstate. The radial wave functions have the same number 
of radial quantum numbers and, again, the lower components agree better. 

Next we examine the other category of Dirac eigenstates with j = I — 1/2 and wave 
functions as in Eq. (3b). Their partners with j' = j + 1 are obtained from Eqs. (8b, c). As 
an example, we consider the realistic relativistic mean field eigenstates Osi/2, lsi/2, 2si/2 
[I = 1, j = 1/2) for 208 Pb [|10j. The Osi/2 eigenstate will have a partner which we de- 
note by [P(0si/2)]d 3 /2, but there is no d 3 / 2 eigenstate at approximately the same energy 
as the Osi/2 eigenstate, so there is no eigenfunction to compare to. On the other hand, 
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the lsi/2 and 2si/ 2 eigenstates are almost degenerate with the 0^3/2 and Ig^/2 eigenstates 
respectively. In Fig. 3 we compare the spatial wave functions of the pseudospin partners 
[P(lsi/ 2 )]d 3 / 2 , [-P(2si/ 2 )]d3/2 with the respective 0d 3 / 2 , ld s / 2 eigenstates. These partners 
agree well with the eigenfunctions in the interior but not on the nuclear surface. In fact, the 
partners do not have the same number of nodes and do not fall off exponentially but inversely 
as the cubic power, r~ 3 , in agreement with the r~( £+2 ) behavior reported in Eq. (8b). 

The Dirac eigenstates with n r = and j = £ — 1/2 are special, because no eigenstates 
exist with the quantum numbers of their partners, h r = and j = £ + 1/2. An example 
is given in Fig. 3a,b for £ = 1, j = 1/2. For heavy nuclei these states with large j are 
the "intruder" states. Before the SU(2) algebra of pseudospin was discovered, these states 
were discarded from the pseudospin scheme. However, that is clearly not a valid procedure 
if pseudospin symmetry is a symmetry of the Dirac Hamiltonian. As another example, we 
show in Fig. 4a,b the radial wavefunction of the [^(0/7/2)] ^9/2 partner of the O/7/2 intruder 
state (£ = 4, j = 7/2). There is no quasi-degenerate h 9 / 2 eigenstate to compare to. The 
upper component has the r" 6 falloff alluded to above. Although both components have zero 
radial quantum number, they do not compare well with the 0h 9 / 2 eigenstate shown in Fig. 
4c, d. In Fig. 4c, d we show also the radial wavefunction of the [P(lf 7 / 2 )}h$/ 2 partner of the 
I/7/2 state {£ = 4, j = 7/2) and compare it to the 0h$/ 2 eigenstate. The upper component 
has again the r~ 6 falloff and therefore does not compare well on the surface. Also the number 
of radial quantum numbers differ. The lower components agree better. 

In summary, we have shown that the radial wave functions of the upper components of 
the j = £ — 1/2 pseudospin partner of the eigenstate with j = £ + 1/2 is similar in shape 
to the j = £ — 1/2 eigenstate but there is a difference in magnitude. However, the n r ^ 
radial wave functions of the upper components of the j = £ + 1/2 pseudospin partner of the 
eigenstate with j = £ — 1/2 is not similar in shape to the j = £+1/2 eigenstate. In fact these 
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wave functions approach r £_1 rather than r £+1 for r small, do not have the same number of 
radial nodes as the eigenstates, and do not fall off exponentially as do the eigenstates, but 
rather fall off as r~( e+2 \ Furthermore, the pseudospin partners of the "intruder" eigenstates, 
h r = 0, also fall off as as r ~C~ +2 \ We have confirmed that the radial wave functions of the 
lower components of the pseudospin partners of eigenstates of the Dirac Hamiltonian for 
j — £±1/2 are very similar to the eigenstates with the same quantum numbers except for 
some differences on the surface. 

This research was supported in part by the United States Department of Energy under 
contract W-7405-ENG-36 and in part by the U.S. -Israel Binational Science Foundation. 
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FIGURES 



Figure 1. a) The upper component [g(r)} and b) the lower component [f(r)] in (Fermi) -3 / 2 
of the [P(0c?3/2)]si/2 partner of the Ods/2 eigenstate compared to the lsi/2 eigenstate. c) The 
upper component and d) the lower component of the [P(ld 3 /2)]si/ 2 partner of the 
eigenstate compared to the 2si/ 2 eigenstate for 208 Pb [HI. 



Figure 2. a) The upper component [g(r)] and b) the lower component [f(r)] in (Fermi) 3 / 2 
of the [-P(0/i 9 / 2 )]/7/2 partner of the 0/i 9 / 2 eigenstate compared to the I/7/2 eigenstate for 



208 



Pb gO 



Figure 3. a) The upper component [g(r)} and b) the lower component [f(r)] in (Fermi) -3 / 2 
of the [P(0si/ 2 )]d 3 /2 partner of the 0si/ 2 eigenstate. c) The upper component and d) the 
lower component of the [-P(lsi/2)]^3/2 partner of the lsi/2 eigenstate compared to the 0^3/2 
eigenstate. e) The upper component and f ) the lower component of the [P(2si/2)]c^3/2 partner 
of the 2sx/ 2 eigenstate compared to the 1^3/2 eigenstate for 208 Pb [[10 



Figure 4. a) The upper component [g(r)} and b) the lower component [f(r)] in (Fermi) -3 / 2 
of the [-P(0/7/ 2 )]^g/2 partner of the 0/ 7 / 2 eigenstate. c) The upper component and d) the 
lower component of the [P(l/ 7 / 2 )]/i9/2 partner of the I/7/2 eigenstate compared to the 0h 9 / 2 
eigenstate for 208 Pb |T0]. 
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